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In this paper, dedicated to the career of Tom Erber, we consider the Casimir interaction between 
weakly coupled bodies at nonzero temperature. For the case of semitransparent bodies, that is, ones 
described by (5-function potentials, we first examine the interaction between an infinite plane and an 
arbitrary curved surface. In weak coupling, such an interaction energy coincides with the exact form 
of the proximity force approximation obtained by summing the interaction between opposite surface 
elements at arbitrary temperature. This result generalizes a theorem proved recently by Decca et 
al. We also obtain exact closed-form results for the Casimir energy at arbitrary temperature for 
weakly-coupled semitransparent spheres. 
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I. INTRODUCTION 



Since the earliest calculations of fluctuation forces between bodies [lj, that is, Casimir or quantum vacuum forces, 
multiple scattering methods have been employed. Rather belatedly, it has been realized that such methods could be 



used to obtain accurate numerical results in many cases [2H5J. These results allow us to transcend the limitations of 
' the proximity force theorem (PFT) [(| 0|, and so make better comparison with experiment, which typically involve 
(~~. , curved surfaces. (For a review of the experimental situation, see Ref. (8|.) 

0^ ■ The multiple scattering formalism, which is in principle exact, dates back at least into the 1950s Particularly 
noteworthy is the seminal work of Balian and Duplantier [ll|. (For more complete references see Ref. [l2| .) This 
technique, which has been brought to a high state of perfection by Emig et al. [5j, has concentrated on numerical 
results for the Casimir forces between conducting and dielectric bodies such as spheres and cylinders. For recent 
impressive numerical results for metals and dielectrics see Refs. [H, [H|]. Our group has noticed that the multiple- 
scattering method can yield exact, closed- form results for bodies that are weakly coupled to the quantum field (I3.fl5j. 
, (That is, we are carrying out first-order perturbation theory in the background potential. For early examples of this 
in the Casimir context, see Ref. [Hj].) This allows an exact assessment of the r ang e of applicability of the PFT. 
The calculations there, however, as those in recent extensions of our methodology have been restricted to scalar 
fields with (5-function potentials, so-called semitransparent bodies. (These are closely related to plasma shell models 
[3.IT3-21].) The technique was recently extended to dielectric bodies [22ll23l |. characterized by a permittivity e. Strong 
coupling would mean a perfect metal, e — > oo, while weak coupling means that e is close to unity. 

In this paper we will extend the weak-coupling formalism to the situation of nonzero temperature. This extension 
is extremely straightforward. We then apply the general formula to the case of an arbitrarily curved semitransparent 
surface above an infinite semitransparent plane. Remarkably, the result coincides with the use of the so-called 
proximity force approximation (PFA), which in its general form is exact in this case for all separations between the 
surfaces and for all temperatures. We also obtain exact closed-form results for the forces between separated spherical 
shells for all temperatures. In the Appendix we discuss exact formulas for arbitrary positive and negative powers of 
the distances between points on two spheres, needed for such calculations. 
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II. MULTIPLE SCATTERING DERIVATION OF VACUUM ENERGY BETWEEN WEAKLY 

COUPLED POTENTIALS 

The quantum vacuum energy for the interaction mediated by a massless scalar field between two nonoverlapping 
potentials Vi(x) and V-^x) is 

E = -~Trln(l-V 1 GiV 2 G 2 ), (2.1) 
in terms of the single potential Green's functions 

Gi = (1 + GoV^Go. (2.2) 

The free Green's function, satisfying 

-d 2 G Q (x-x') = 5(x-x'), (2.3) 

has the explicit form 

G (x-x')= f ^ (r-r , , W )e-»('- i '), (2.4) 



where the time-Fourier transform is 



e -|Clk-r'| 

S°( r - r ''«C) = 4d737T (2 - 5) 



where we have performed the Euclidean rotation uj — > i(. 
For weak potentials, the energy (|2.1|) simplifies dramatically: 



E « ^rV lGo V 2 G = * [ (dr)( d r')*mn. (2.6) 



— 647r 3 

At finite temperature the integral over imaginary frequency becomes the Matsubara sum: 

ja „— 2|f I |r— r'| 00 — inT\m\ |r— r'| rp 

7vT^ T Y ~\ n2— = 1 ^coth27rT|r-r'|, (2.7) 



— oo 



m— — oo 



so the interaction energy becomes 

T fr^\rj„>w ^t/ ^coth27rT|r - r' 

327T 



^ = / (dr)(dr')Vi(r)F 2 (r') 1^ (2.. 



which evidently reduces to Eq. (|2.6p for T = 



III. PARALLEL PLATES 

For parallel, semitransparent plates, separated by a distance a, where the potentials are 

Ki(r) = Ai«(z), V 2 (r) = X 2 5(z-a), (3.1) 
the integrals in Eq. (|2.6|) are readily carried out, with the resulting energy per unit area A, £ = E J A: 

327r z a 

This well-known result holds even if one of the plates has a finite area A. At finite temperature the result is 

AiA 2 T f°° dx , , , 

ct = / — cothx. (3.3) 

16tt J 27rTa sc 

The energy is ambiguous because it depends on the arbitrarily chosen upper limit. However, it corresponds to a 
well-defined pressure between the plates, 

P T = ~&r = -%^coth2^Ta. (3.4) 
da 167ra 



3 



IV. INTERACTION BETWEEN AN INFINITE PLANE AND AN ARBITRARILY CURVED SURFACE: 

PFA 

Now consider the interaction between a semitransparent plane, described by the potential 

Vx(r) = \ 1 6(z), (4.1) 

and an arbitrary curved surface S, which does not intersect the plane z = 0, which corresponds to the potential 

V 2 (r) = \ 2 S(z-s(x,y)), (4.2) 

where z = s(x, y) is the equation of the surface. Then, from Eq. (|2.8[) it is immediate that the energy is (the upper 
limit of the x integration is again physically irrelevant) 



E T = -h^L[ dS { dx c J^ i (4.3) 

J2ttTz(S) X 



16tt 



where the area integral is over the curved surface. This is precisely what one means by the proximity force approxi- 
mation, where one sums energies between adjacent elements treated as parallel plates: 

£ PFA = / dSE\\(z(S)), (4.4) 

in view of Eq. (|3.3p . This is in fact just the theorem proved by Decca et al. 124], who were considering gravitational 
and Yukawa type forces, but we see it applies to any central force. 

For example, the above, exact formula for weakly-coupled semitransparent surfaces says that the force on such a 
sphere, of radius a, the center of which is a distance Z above a semitransparent plane is 

dE T X^aT r^nz+a) du 
F T = -—- = / — cothu. (4.5) 

The zero-temperature limit of this is 

87T Z z — a z 

which may be alternatively derived from the zero-temperature energy 

_ AiA 2 a 2 f 1 dcos9 AiA 2 a Z + a 

16tt J_ X Z + aco*6~ 16tt n Z - a ' ' ' 

again, the exact PFA result. 



V. INTERACTION BETWEEN TWO SEMITRANSPARENT SPHERES AT NONZERO 

TEMPERATURE 

Consider now two spheres, of radius a and b, respectively, with a distance between their centers R > a + b. In terms 
of local coordinates with origins at the centers of the two spheres, the semitransparent potentials are 

Vi = \i5(r-a), V 2 = X 2 S(r'-b), (5.1) 

and let us further suppose that R lies along the z axis of both coordinate systems. Then the squared distance between 
points on the spheres is 

|r - r'| 2 = R 2 + a 2 + b 2 - 2abcos~f - 2R(acos9 - bcos9'), (5.2) 

in terms of polar angles in the two spheres, where the cosine of the angle between the two radial vectors locating the 
points is 

cos 7 = cos 9 cos 9' + sin 9 sin 9' cos (</> — <//). (5-3) 
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We insert this into the expression for the energy (|2.8[) . obtaining 

coth27rT|r- r'| 



E 



XlX2T a 2 b 2 I dndn' 



32tt 2 

It seems difficult to proceed in general, but we can work out a low temperature expansion using 

00 o2n l > 111 
1 -°2n 2n-l _ 1 , 1 1 



COthy = £ ^-^V™- 1 = ± + ±y - + 

' [2n)\ v 3 45 

n=0 v ; y 

which will give rise to an expansion of the form 

E T = E Q + T 2 E 2 + T 4 E A + . . 
The zero temperature result was worked out, by inspection, in Ref. (T^ |: 



(M <0 



£ 



AiA 2 a& , 1 - (a - 6) 2 /i? 2 
■ In ■ 



WttR l + {a + b) 2 /R 2 ' 



(5.4) 



(5.5) 



(5.6) 



(5.7) 



The T 2 term is trivial because it is evaluated by Newton's theorem that a Coulomb potential exterior to a spherically 
symmetric charge distribution is as though the charge were concentrated at the center: 



E 2 = - 



\1\2n a 2 b 2 
~3 R~ 



The T 2n term, n > 1 however, is slightly nontrivial: 



„ AiA 2 2u2 {^) 2n B 2n /\ , ,, 



We may evaluate the integrals by expanding in powers of a = a/R and b = b/ R: 

1 



dn<m'\v-r'\ = (4n) 2 R 
dSldCl' \r-r'\ 3 = {Att) 2 R 3 
dn dQ'Ir-r'l 5 



l + -(a 2 +b 2 ) 



(4nfR 5 
1 



1 + 2(a 2 + b 2 ) + ia 4 + -d 2 b 2 + -b 4 
5 3 5 

1 + 5(a 2 + b 2 ) + 3a 4 + 10a 2 6 2 + 3& 4 



-a 6 + a 2 b 2 (a 2 + b 2 ) + -b 6 
7 7 



(5.8) 
(5.9) 

(5.10a) 
(5.10b) 

(5.10c) 



and so on. The reason these are polynomials is evident when one considers the multipole expansion of the Coulomb 
potential — See, for example, Chap. 22 of Ref. |25j. For general formulas for such moments, see the Appendix. 
By computing further terms in the sequence of polynomials, we are able to recognize the pattern: 



• „ n+l p 



(47T) 



rXp-Q)j 



=0 9=0 



where 



A(n,p,q) 



(2n + 2)\ 



(2n -2p + 2)!(2p -2q+ l)\(2q + 1)! ' 
When this is inserted into the low temperature expansion, we can remarkably sum the series: 



(5.H) 



(5.12) 



En 



AiAaofef I -{a- bf / R 2 
1671? \ 11 1 - (a + b) 2 /R 2 



f(2irT(R + a + b)) + f{2i:T{R - a - b)) - f(2nT(R - a + b)) - f(2irT(R + a-b))\, (5.13) 



5 



where / is 



M = E 



..2n 



^ 2n(2n- l)(2n)r ' 
which is obtained from the second antiderivative of the hyperbolic cotangent: 

i_/( !/ ) = cothy-i /(0)=/'(0) = 0. 



dy 



y 



(5.14) 



(5.15) 



Although the power series expansion (|5. 14)) is valid only for sufficiently low temperatures 2T(R + a + b) < 1, the 
solution of the differential equation is valid for all values of T. 

For sufficiently high temperatures we can replace the hyperbolic cotangent in the differential equation by 1, and 
then 



f(y) ~y\ny + lny + Ay + B, y > 1, 



(5.16) 



where A and B are integration constants that do not contribute to Eq. (|5.13|) . When this asymptotic solution is 
inserted into Eq. (I5.13[) the zero temperature logarithm cancels out, and we are left with 



Ex 



\1\2ab 



T 



In 



R 2 -(a + bf 
R 2 -(a- by 



a (R + b) 2 - a 2 
R n (R - b) 2 - a 2 



This result may be derived directly from the high-temperature form 



En 



X 1 X 2 Ta 2 b 2 
32^ 



dn dn'- 



1 



b_ (fi + a) 2 - b 2 
R n (R - a) 2 - b 2 



T 



T 



00. 



(5.17) 



(5.18) 



This again may be worked out by expanding in powers of the radii of the spheres. Computing the first several terms 
reveals the pattern: 



dn dfl'-, 



1 



r'|2 



R 2 



E 



1 



(n + l)(2n + l)f r 



1) ^-C 2 



1 f 2n + 2 
2m + 1 



i 2 < 



n — 7n) Lzm 



(5.19) 



This sum is almost identical to that found for spheres at zero temperature, as seen in Eq. (6.15) of Ref. pJl, which 
led to Eq. (|5.7p . except for the appearance of l/(2n+ 1) here. Therefore, the former series must be obtained from the 
present series by differentiation. Denoting the double sum in Eq. (|5.19[) by S, it must be true that 



R z 



d S R 2 



In 



1 - (a + b) 2 /R< 



(5.20) 



dRR 4ab*"\l-(a-b) 2 /R 2 , 

where S is R 2 /4ab times the square-bracketed quantity in Eq. (|5.17[) . This equality is, in fact, easily verified. Sec the 
Appendix for the generalization of this result. 

We compare the general form [obtained by numerically integrating Eq. (|5 . 1 5[) ] and the high-temperature limiting 
form (|5TT7|> in Fig. [TJ 



VI. CONCLUSIONS 



We have shown that exact results may be found in weak coupling for the quantum vacuum forces between nontrivial 
bodies not only at zero temperature, but at finite temperature. We have shown that the exact form of the proximity 
force approximation holds exactly for all temperatures for the force between an infinite plane surface and an arbitrarily 
curved one. We have also computed the force between two semitransparent spheres at arbitrary temperatures, and 
obtain remarkably simple, closed-form expressions. The PEA equivalence evidently will hold for tenuous dielectric 
bodies in electromagnetism, and closed-form finite temperature results may be easily obtained between dielectric 
bodies as well. 
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FIG. 1: Comparison between the general and high temperature forms of the energy, as a function of RT. Energies are shown 
for a — b = R/4. The high temperature result is linear in T. Also shown is the power series expansion [Eq. (|5.14|) truncated at 
200 terms], which diverges in this case at RT =1/3. Plotted is e = — 167r_RE/(AiA2a 2 ). 



Appendix A: Mean powers of distances between points on spheres 

In Sec. [V] we used exact evaluations of mean distances, defined by 

dfldn'\r -r'\ p = (AirflFPpia,!)), 



(Al) 

where R is the distance between the centers of the two nonoverlapping spheres, of radii a and b, respectively. Here 
a = a/R and b — b/R, and P p (a, b) can in general be represented by the infinite series 



^(^)=g (2 „ 2 +2)! r( r ( "_;-T) 1) g - (a ' S) - 

Here the homogeneous polynomials Q n are 



= 1, 

= 2(a 2 + 6 2 ), 

= 3a 4 + 10a 2 ,? + 36 4 , 

= 4a 6 + 28a 4 6 2 + 28a 2 6 4 + 4b 6 , 



or in general, 



n 



2n + 2 
2m + 1 



^2(n — ra)yi7n 



We can easily see the following recursion relation holds: 



P p -^a,b) = ^--^R^P p (a 1 b), 

since Q n is homogeneous in R of degree — 2n. 

For p a non- negative integer, P p is a polynomial of degree 2|~p/2], and we can immediately find 



(A2) 



(A3a) 
(A3b) 
(A3c) 
(A3d) 



(A4) 



(A5) 



P p (a,b) 



1 



4a6(p + 2)(p + 3) 



(1 + a + b)' p+3 + (1 - 6 - b) p+3 - (1 - a + b) p+3 - (1 + a - 



p = 0,l,2, 



(A6) 
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For p a negative integer, we have 



P-i 



1. 



(A7a) 



P-2 




(1 + 6) 2 

(1-6)2 



— a 



+ 6 In 



(l + a) 2 -6 2 



(1 - a) 2 - 6 2 



(A7b) 



P-3 



1 , l-(a + 6) 2 

-T7? ln : — n — 



4a6 1 - (a - 6) 2 



(A7c) 



P-4 



1 



(A7d) 



[l-(a + 6) 2 ][l-(a-6) 2 ]' 



and further expressions can be obtained by use of Eq. (|A5|) . 
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